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Abstract: We compute conductivities of strongly-interacting and non-uniform charge den¬ 
sities dual to inhomogeneous anti-de Sitter-black hole spacetimes. Backreacting bulk scalars 
with periodic boundary profiles, we construct generalizations of Reissner-Nordstrom-AdS that 
interpolate between those used in two previous studies — one that reports power-law scal¬ 
ing for the boundary optical conductivity and one that does not. We find no evidence for 
power-law scaling of the conductivity, thereby corroborating the previous negative result that 
gravitational crystals are insufficient to generate the power-law mid-infrared conductivity 
observed in cuprate superconductors. 

Keywords: Gauge/gravity duality, holography, condensed matter physics (AdS/CMT) 
ArXiv ePrint: 1506.06769 

1 Guggenheim Fellow 




Power laws are ubiquitous to critical phenomena as they are the fingerprint of scale invari¬ 
ant correlations. To no surprise then, quantum criticality is the most commonly proffered 
explanation[l-3] for the power-law scaling 

|<r(o;)| ~ cu -2 / 3 (0.1) 

observed in the mid-infrared frequency range of the optical conductivity[4-7] for optimally 
doped copper-oxide superconductors. However, we currently have no microscopic under¬ 
standing of how quantum criticality emerges from the strong correlations that mediate the 
normal state near optimal doping. Such an understanding requires precise knowledge of the 
low-energy degrees of freedom in the strongly coupled regime. Ascertaining these degrees of 
freedom has proven difficult because the integral of the optical conductivity [8, 9] up to the 
optical gap exceeds the number of doped holes. Consequently, no one-to-one mapping[10] 
exists between the number of doped holes and the actual propagating degrees of freedom at 
low energy. 

Given the seemingly unbridgeable divide between the ultraviolet (UV) and infrared (IR) 
physics in the cuprates, it is desirable to address the problem with a method which is based on 
conserved currents rather than on a traditional particle description. The gauge/gravity dual¬ 
ity offers such an alternative, in which a strongly-coupled, conformal, d-dimensional quantum 
field theory lives at the boundary of a (d + l)-dimensional gravity theory. The workhorse 
in condensed matter applications is Reissner-Nordstrom-AdS (RN-AdS), which provides the 
simplest bulk description of matter at finite charge density. RN-AdS is often treated as a 
static background to be augmented by a set of probe fields which source additional operators 
at the boundary. Horowitz, Santos, and Tong (HST)[11] improved upon these treatments 
by solving the completely coupled equations of motion for gravity, electromagnetism, and 
a neutral scalar field. Further, they fixed the boundary source provided by the scalar to be 
Aq cos kx, thereby imprinting a periodic distortion upon the bulk geometry and the boundary 
charge density. Taking this inhomogeneous bulk as a new static background, the boundary 
conductivity was then obtained by solving equations of motion for propagating fluctuations 
of the bulk fields. The key finding of reference [11] is that the conductivity obtains — in 
addition to a Drude peak — a mid-infrared power law identical to that seen in the cuprates, 
eq. (0.1), plus an additive constant. HST obtained identical power laws for an ionic lattice — 
charged matter with a spatially periodic chemical potential, lirn^o At = fj,(x) — in a slightly 
different geometry[12]. These conclusions are indeed startling and imply that a gravitational 
crystal encodes the optical conductivity of the cuprates. 

However, the key claim that Einsteinian gravity, a Maxwell held and an inhomogeneous 
charge density encode the mid-infrared conductivity of the cuprates has been called into 
question[13, 14]. In reference [13] Donos and Gauntlett (DG) studied a model inspired by 
the Q-ball potential of Coleman[15] which has the added simplification that only ordinary 
rather than partial differential equations need be solved to obtain the conductivity. Their 
scalar held, of the form <h(z, x) = cj)(z)e lkx , leads to a uniform charge density as they chose a 
potential of the form E(|<h| 2 ). Added differences with the work of HST is the use of a scalar 
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mass of m 2 = —3/2 1 and the radial gauge as opposed to the de Donder gauge and Lorenz 
gauges. In addition, rather than using a log-log plot to discern the presence of a power law, 
they plotted 

\&\" 

a = 1 + uj- —(0-2) 

\a\ 

A value of a = —2/3 would correspond to the mid-infrared power-law conductivity of the 
cuprates. None was found. In fact, a was found to vary fairly widely as a function of 
frequency precisely where the power law was reported by HST. Similar results from a slightly 
different model have also been reported in reference [14]. 

Given the substantial differences between the HST and DG constructions, it is instructive 
to resolve this problem through a new construction that interpolates between the two models. 
We achieve this by introducing two neutral scalars, subject to periodic potentials with tunable 
phase angles, to four-dimensional RN-AdS. Although the HST and DG setups differ drastically 
in terms of the boundary charge density, we find that for identical system parameters used by 
HST, a deviates significantly from the constant needed to reproduce the mid-infrared power 
law of the cuprate conductivity. Consequently, we conclude that Einsteinian gravitational 
crystals are insufficient to explain the power-law scaling in the cuprates. 

We consider the Einstein-Maxwell action 

S =usk/ A ^( R - 2 A -H' (a3) 

where A = —6/L 2 is the cosmological constant, L is the AdS length, F = dA is the field 
strength of a Maxwell field, R is the Ricci scalar, g is the metric determinant, and Gjv is the 
Newton constant. To this we append an action for two neutral scalar fields, 


S* = 


1 


167rG/v 


d 4 x^Yl [2 (V<D 4 ) 2 + 4R (<!>,,) 




(0.4) 


The asymptotic expansion for the scalar fields is 


} + z A <f >{ 2]1 + 

^(x) = Aicos ( kix - ^ ) , 


4>^(x) = A 2 cos ^k 2 x + 0 , 


(0.5) 


1 Extremal RN-AdS4 has an emergent AdSi x R 2 geometry near the horizon, which hosts a quantum dual 
with a Breitenlohner-Freedman (BF) bound higher than that of the boundary theory. Donos et al. argue 

that the HST results may be unstable because their chosen scalar field mass violates the BF bound in the 
near-horizon CFT. 
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where A = 3/2 + \/9/4 + m 2 . 

Obtaining the boundary conductivity requires two distinct calculations. First, solve the 
full equations of motion for a static background. Second, kick the system with a time- 
dependent electric field and find the current response. As all of this is standard (see references 
[11, 12, 16, 17]), we will only highlight the key features. To assist in solving the equations of 
motion, 


Eab = Rab - Rgab ~ Tab = 0 , ( 0 . 6 ) 

V a F ab = o, (0.7) 

□$i - v 1 ($i) = 0, (0.8) 

we solved instead the Einstein-DeTurck equation 

= (0.9) 

f“ = s“‘(r£,(«)-r£,(5)). (o.io) 


Here T a j, = T a b — {g a b/ 2 ) T is the trace-reversed energy-momentum tensor and g is a reference 
metric. To directly compare with the HST results, we use the mass value m 2 = —2. The 
Einstein-DeTurck equation simultaneously encodes E a b = 0 and the harmonic gauge £“ = 0 
due to E a b obeying a continuity equation. Because the DeTurck term breaks gauge invariance, 
metric components can be turned on to compensate any terms from T a ^. We take the ansatz 


ds 2 = —7T 


2 r 


2 2 L 


- (1 - z)P(z)Q a dt 2 + 


Qzzdz 


(1 - z)P(z) 
pQxx(dx T z Qzxdzj -{- Qyydy , 

P{z) = 1 + z + z 2 - ^z 3 , 

A = (1 — z)at(z, x)dt, <1 \ = z(f>i(z,x). 


( 0 . 11 ) 


Here all coordinates have been rescaled to be dimensionless and the radial coordinate z is 
parameterized to extend from [0,1]. The conformal boundary exists at z = 0 and an event 
horizon exists at z = 1. The temperature is given by T = P(1)/(AttL). This ansatz clearly 
reduces to Reissner-Nordstrom-AdS (RN-AdS) when Q tt = Qzz = Qxx = Qyy = 1 ,Q Z x = 0, 
at = gi = n and <fii = 0, which is used as our reference metric g in (0.10). 

If we set A\ = A 2 = Ao and k\ = k 2 = k , then the parameter 9 tunes this model between 
the HST and Q-lattice models. 9 = 0 corresponds to the HST lattice and 9 = ir/2 yields the 
Q-lattice. Consequently, regardless of the origin of the power law, the model considered here 
should be able to unearth the source. 

We have eight equations of motion to solve. At the conformal boundary, z = 0, we 
imposed Dirichlet boundary conditions. These conditions will produce the desired AdS ge¬ 
ometry, fix the chemical potential, and fix the scalar lattice. The remaining conditions will 
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6»=tt/2 (Q) 
(2=7t/4 

6>=0 (HST) 



Figure 1: (Left) Plot of Q zx for Aq = 0.75/\/2, k = 1, 0 = 0, /r = 1.4, and T / p, = 0.115. 
For this plot max(|£ a |) < 10“ 12 . (Right) Plot of the charge density p = lim^o \J~9F tz for 
same parameters, and with 6 = 7t/4, ir/2. 


be so-called regularity conditions. We demand that all the functions can be expanded at the 
horizon z = 1 as 


Qab(z,x ) = Q ( a 0 b \x) + (1 - z)Q ( a l\x) H -, 

a t (z,x ) = a[°\x) + (1 — z)a[ 1 \x) H-, 

(,k(z,x ) = </>- 0) (x) + (1 - z)4>f > {x) H-. (0-12) 

We then plug these expansions into various equations. Two further conditions can be specified 
by setting £ z = £ x = 0 and the remaining six will be set by demanding that the equations of 
motion be satisfied to lowest order in the expansion. This will fix relationships among the (0) 
and (1) coefficients, giving Robin type boundary conditions at z = 1. The numerical solution 
was obtained through the Newton-Raphson method using a sixth-order finite difference grid 
for the z-direction and a Fourier grid for the x-direction. Details of such methods can be 
found in reference [18]. 

In figure 1 we have a sample plot of the off-diagonal component Q zx ■ It is clear that the 
periodicity of the scalar field is imprinted on the background (in multiples of 2k due to the 
quadratic form of the scalar field in the energy-momentum tensor). Figure 1 also contains 
the boundary charge density for different values of 6, distinguishing the periodic imprint of 
the HST lattice versus the uniformity of the Q-lattice. Figure 1 was computed in MATLAB 
at double precision with (300,45) grid points on the (z, x) plane, which is a typical grid size 
for the calculations in this paper. 

The second step involves a time-dependent perturbation around the static solution which 
enables a calculation of the conductivity. We follow the methods of reference [17]. Denoting 
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Figure 2: Low frequenecy plots of real and imaginary parts of the conductivity for various 
parameters. In each plot Ao/k = 3/(4\/2), fi = 1.4 and T/^l = 0.115. 
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Figure 3: Plots of the power function and the 
parameters as figure 2. The dotted line marks 
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argument of the conductivity for the same 
—2/3. Refer to figure 2 for legend. 


background fields with bars, we write 


Sab — Sab ~b habi A a — A a T b a , T (0.13) 

where the barred quantities are the static background and the extra pieces h a b , b a , and iji 
are fluctuations. The leading term in b x sources an electric field at the boundary and the 
subleading term contains the current response. The Kubo formula can then be used to obtain 
the conductivity. 

Using this method, we computed the conductivity as a function of the interpolating 
parameter 9. Shown in figure 2 are the real and the imaginary parts of the conductivity for 
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Figure 4: This plot was generated using A\ = 0.75, k\ = 1, = 2, 0 = 0, fi = 1.4 and 

T /ju = 0.115. The amplitude A 2 of the higher harmonic lattice was adjusted. 


two different values of k and three values of 0. For each choice of parameters the low frequency 
conductivity obeys the Drude form. Note, even though 9 = 7r/2 corresponds to a uniform 
charge density of DG, the conductivities are almost identical to those of HST (0 = 0). 

Figure 3 is the key test for the existence of the power-law conductivity. In the left panel 
we show a couple of plots of a using the parameter choices of reference [11]. The dotted 
line indicates a power law of —2/3. As is evident, regardless of the value of 9, no discernible 
power law exists even as k is varied. Also of note is the fact that the DG (0 = vr/2) and HST 
(0 = 0) models yield almost identical numerical results for the conductivity. The right panel 
of figure 3 presents the phase angle which also deviates from 60°. 

As a final check we also varied the amplitude A 2 for a dichromatic lattice in figure 4, as 
this could introduce a mix of higher harmonics. In this case as well, no evident power law 
exists. In fact, for any range of parameters including temperature, no power law was found. 

Since we have introduced a model that is capable of interpolating between DG and HST 
and we find no power law in either case, we conclude that gravitational crystals, although 
adequate in describing Drude response, do not encode the power-law optical conductivity of 
the cuprates. 

Regarding the origin of the power-law optical conductivity, the only study[19] to date 
that has successfully reproduced the u;~ 2 / 3 scaling relies on excitations which exist on all 
energy scales — namely scale-invariant matter or unparticles. Given that the radial direction 
in AdS represents the running of coupling constants, in principle it contains the correct 
ingredients to capture unparticle excitations. Hence, we anticipate that some construction 
using gauge/gravity duality, other than the one presented here, should be able to reproduce 
the power law. 


- 7 - 













Acknowledgments 


We thank Aristos Donos and Jorge Santos for their generosity in assisting us with the 
numerics. P. Phillips thanks the NSF DMR- 1104909 and the J. S. Guggenheim Foundation 
for financial support. 

References 

[1] P. W. Anderson, Infrared conductivity of cuprate metals: Detailed fit using Luttinger-liquid 
theory, Phys. Rev. B 55 (May, 1997) 11785-11788. 

[2] M. Ogata and P. W. Anderson, Transport properties in the Tomonaga-Luttinger liquid, Phys. 
Rev. Lett. 70 (May, 1993) 3087-3090. 

[3] E. Abrahams and C. M. Varma, What angle-resolved photoemission experiments tell about the 
microscopic theory for high-temperature superconductors, Proceedings of the National Academy 
of Sciences 97 (2000), no. 11 5714-5716. 

[4] D. v. d. Marel, H. J. A. Molegraaf, J. Zaanen, Z. Nussinov, F. Carbone, A. Damascelli, 

H. Eisaki, M. Greven, P. H. Kes, and M. Li, Quantum critical behaviour in a high-Tc 
superconductor, Nature 425 (2003), no. 6955 271-274. 

[5] A. El Azrak, R. Nahoum, N. Bontemps, M. Guilloux-Viry, C. Thivet, A. Perrin, S. Labdi, Z. Z. 
Li, and H. Raffy, Infrared properties ofYBa 2 Cu 30 i and Bi 2 Sr 2 Ca n -\Cu n 02 n +A thin films, 
Phys. Rev. B 49 (Apr, 1994) 9846-9856. 

[6] J. Hwang, T. Timusk, and G. D. Gu, Doping dependent optical properties of Bi 2 Sr 2 CaCu 20 s+, 
Journal of Physics Condensed Matter 19 (Mar., 2007) 125208, [cond-mat/0607653]. 

[7] Z. Schlesinger, R. T. Collins, F. Holtzberg, C. Feild, S. H. Blanton, U. Welp, G. W. Crabtree, 

Y. Fang, and J. Z. Liu, Superconducting energy gap and normal-state conductivity of a 
single-domain YBa 2 CusOT crystal, Phys. Rev. Lett. 65 (Aug, 1990) 801-804. 

[8] S. Uchida, T. Ido, H. Takagi, T. Arima, Y. Tokura, and S. Tajirna, Optical spectra of 
La 2 - x Sr x CuO 4 : Effect of carrier doping on the electronic structure of the CUO 2 plane , Phys. 
Rev. B 43 (Apr, 1991) 7942- 7954. 

[9] Y. S. Lee, K. Segawa, Z. Q. Li, W. J. Padilla, M. Dumm, S. V. Dordevic, C. C. Homes, 

Y. Ando, and D. N. Basov, Electrodynamics of the nodal metal state in weakly doped high- T c 
cuprates, Phys. Rev. B 72 (Aug., 2005) 054529, [cond-mat/0509167]. 

[10] P. Phillips, Colloquium: Identifying the propagating charge modes in doped Mott insulators, Rev. 
Mod. Phys. 82 (May, 2010) 1719-1742. 

[11] G. T. Horowitz, J. E. Santos, and D. Tong, Optical conductivity with holographic lattices, 
Journal of High Energy Physics 7 (July, 2012) 168, [arXiv: 1204.0519]. 

[12] G. T. Horowitz, J. E. Santos, and D. Tong, Further evidence for lattice-induced scaling, Journal 
of High Energy Physics 11 (Nov., 2012) 102, [arXiv: 1209.1098]. 

[13] A. Donos and J. P. Gauntlett, Holographic Q-lattices, Journal of High Energy Physics 4 (Apr., 
2014) 40, [arXiv: 1311.3292], 


[14] M. Rangamani, M. Rozali, and D. Smyth, Spatial modulation and conductivities in effective 
holographic theories, Journal of High Energy Physics 7 (July, 2015) 24, [arXiv: 1505.05171]. 

[15] S. Coleman, Q-balls, Nuclear Physics B 262 (1985), no. 2 263 - 283. 

[16] M. Headrick, S. Kitchen, and T. Wiseman, A new approach to static numerical relativity and its 
application to Kaluza-Klein black holes, Classical and Quantum Gravity 27 (2010), no. 3 
035002. 

[17] A. Donos and J. P. Gauntlett, The thermoelectric properties of inhomogeneous holographic 
lattices, Journal of High Energy Physics 1 (Jan., 2015) 35, [arXiv: 1409.6875]. 

[18] L. Trefethen, Spectral Methods in MATLAB. Society for Industrial and Applied Mathematics, 

2000 . 

[19] K. Limtragool and P. Phillips, Power-law Optical Conductivity in the Cuprates from Unparticle 
Stuff, ArXiv e-prints (June, 2015) [arXiv: 1506.00649]. 


- 9 - 


